Nonperturbative study of generalized ladder graphs in a (p 2 x theory 
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One of the important issues in the study of a composite 
hadronic system at higher energies, is the search for prac- 
tical and reliable schemes to describe its relativistic dy- 
namics. Our knowledge about the relativistic two-body 
bound state problem in field theory is almost exclusively 
based on the application of the ladder approximation to 
the Bethe-Salpeter equation (BSE) Unfortunately, 
the general applicability of the ladder theory can be ques- 
tioned on physical grounds. In particular, the so-called 
one-body limit does not lead to the Klein-Gordon equa- 
tion as it ought to. Moreover, gauge invariance can not be 
satisfied within this approximation. In order to recover 
these properties, at least the set of all crossed ladder con- 
tributions is needed additionally So far, however, 
the study of the two-body Green function beyond the lad- 
der theory has been considered not feasible in practice. 
With this situation in mind, several quasi-potential equa- 
tions (QPEs) have been proposed and studied as possi- 
ble candidates for an effective theory. Both the ladder 
BSE as well as various QPEs have been used in numer- 
ous studies throughout a wide range of systems, including 
mesons , small nuclei JlO|— |l2| , few-electron atoms jlj] 
and positronium 

In constructing the QPEs, one usually chooses the ap- 
proximations leading to them such that the above men- 
tioned problems are, at least partially, solved. However, 
due to our ignorance of the behavior of the full BSE 
solutions, it is presently unclear which of the, possibly 
infinite number of, QPEs provides the best effective de- 
scription. In this connection it is clearly of interest to 
have actual solutions available for cases where a larger 
class of graphs than the ladder series is included in the 
BSE and that do not suffer from the difficulties inherent 
to the latter approximation. Such solutions may serve as 
a testing ground for the various QPE descriptions. Here 
we present results for the case where in addition also the 
complete set of all irreducible crossed-ladder graphs is 
included in the kernel of the BSE, being the minimal set 
that is free from the above problems. Self-energy and 
vertex corrections are not taken into account. The inclu- 
sion of these contributions are expected not to lead to 
qualitatively different predictions [ 13) . 

In this letter the bound states formed by two scalar 
particles <p with mass m interacting through the ex- 
change of a third scalar particle x with mass /x, are 
determined using the Feynman-Schwinger representation 
(FSR) @Jlf-|lf]. Starting from the Euclidean action for 
the above ip 2 x theory 



S 



A 4 x 



(d x >p) 2 + m 2 ip 2 



+ \{d xX ) 2 + Wx 2 +9^X , (1) 



we may reconstruct the bound state of two (/3-particlcs 
with the set of one-meson-exchange and all irreducible 
crossed-ladder graphs as driving force, by explicitly inte- 
grating out the fields in the two-body Green function G. 



Details of this procedure can be found in Ref. [Q. Ac- 
cording to (J] , the FSR offers a closed expression for the 
'quenched' G (i.e., neglecting the possible occurrence of 
vacuum fluctuation (pip-loops) in terms of path integrals 
over the particle trajectories z and z of the two ip parti- 
cles. Neglecting also the contributions corresponding to 
the self-energy and vertex corrections, it has the form 



G 



As 



As / (Vz) xy / (Vz) s 



exp (-K[z, s] - K[z, s] + V[z, z, s, s}) , (2) 



(3) 



where K and V are given by 
1 f 1 

K[z,s}=m 2 s + — J^At z 2 x {t), 

V[z,z,s,s] = g 2 ss [ At [ Af A (z(t) - z(f)) . (4) 
Jo Jo 

Our main objective here will be to compare the predic- 
tions obtained from Eqs. to those from the ladder 
BSE and various QPEs. 

The functional integrations are over all possible paths, 
subject to the boundary conditions z(0) = x, z(l) = y 
and similarly for z. In (3+1) dimensions the free two- 
point function A(a;) is given by 



A(x) 



4tt 2 | 



;Ki(fi\x\). 



(•5) 



In |jj it was shown that for unequal masses, Eq. (||) satis- 
fies the correct one-body limit. In addition, it was proven 
that, combined with Eq. (Q), it effectively sums up all 
ladder and, due to the absence of any ordering in the 
interaction kernel, also all crossed-ladder contributions 
to G. Each graph of this set is UV-finite, so that no 
short-distance regularization is required. 

The bound state spectrum can be determined by 
studying the behavior of G with respect to variations 
of its initial points (x, x) and final points (y,y). Con- 
sidering, in particular, large timelike separations T = 
h(y4 + §4 — X4 — X4), we infer from the spectral decom- 
position 



G = ^c n exp (-m„T) 



Co exp (— toqT) . 



(6) 



that, asymptotically the Green function is dominated by 
the ground state contribution. 

Notice that the path integrals in Eq. (||) are quantum 
mechanical ones. This amounts to a considerable reduc- 
tion in number of degrees of freedom as compared to, 
for example, putting the field action (Q) on a discrete 
4-dimensional lattice. As a result accurate calculations 
can be carried out with this approach also for very large 
times T. 

Let us now briefly discuss the traditional Bethe- 
Salpeter approach fyfl] to the two-body bound state 
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problem. In the ladder approximation the wave func- 
tion W in momentum space obeys the following integral 
equation 



(27t)< 



dVn?-?')*(?'), (7) 



where q is the relative momentum between the two ip 
particles. After a Wick-rotation, the free two-body prop- 
agator S and the bare interaction V assume the following 
form in the CM-frame 



S(q) = 



V(q - q') = g 



1 



q 2 + Qj 2 + m 2 - js) + soj 2 
1 

(q - q') 2 + (w - Lo') 2 + n 2 ' 



(8) 



(9) 



with the relative momentum q — (q, u>). In the bound 
state region Eq. (Q) only supports solutions for values of 
the invariant energy y/s, that correspond to bound states. 

Since for unequal masses Eq. (^) in the ladder ap- 
proximation does not possess the correct one-body limit, 
several modifications to it have been proposed. Gener- 
ally, they reduce the description from a 4-dimensional 
to a 3-dimensional one by making an ansatz for one of 
the functions involved. This ansatz is chosen such that 
the resulting quasi-potential equation does possess the 
correct one-body limit. Here we study three particular 
examples: the BSLT-equation |lS| ], the equal-time (ET) 
equation j|JT^,|ll|,^] and the Gross equation |5|,[l^] , which 
have been widely used in the literature. 

For the BSLT equation one assumes that the pole 
structure of the two-body propagator can be approxi- 
mated via a dispersion relation. Similar to the BSLT 
case, in the ET prescription the interaction is usually 
supposed to be independent of the relative time, i.e., 
also neglecting retardation effects. An additional term 
is supplied in order to include some of the crossed-box 
contributions. In doing so, the correct one-body limit is 
obtained in this approach. Finally, in the Gross formal- 
ism one puts one of the two particles on its mass-shell by 
hand. These procedures lead to the following forms of 
S(q) 



Sqpe(<7) B = T 2-7T S(u>)- 



2 _i_ m2 n 2 4- m 2 — i.s ' 



\/q 2 + m 2 q 



= 2ir o[u) 



y/q 2 + m 2 q 2 + m 2 - |s 

s 



x 12- 



4(q 2 + m 2 



(10) 



(11) 



Gross 



r-F- (12) 



4 v / s\/q 2 + m2 Vq 2 + m " ~ f V s 

For all cases the delta- function allows for the elimination 
of the relative energy variable to from the description. 



The ladder BSE and 3-dimensional QPEs were solved 
by performing a standard partial wave decomposition, 
thereby factorizing the angular variables. 

The FSR solutions were obtained by discretizing the 
functional integrals, according to 



(Dz) xy 



N 



2N N-l 



n 

i=l 



d 4 ; 



(13) 



The normalization in Eq. (|T^) was chosen such that, 
when expanded in the coupling g 2 , the Green function 
correctly reproduces the Feynman perturbation series. In 
terms of the discretized variables the functionals K and 
V assume the following form 



K[z,s] 
V[z,z,s,a] 



2 ^V^/ 

m s + — 2_^{z t - , 



N 



4.s 



(14) 



2 - N 

g ss s-^ n 



Zj-Zj-i)) (I 5 ) 



The discretized boundary conditions are zq = x, zjq = y 
and similarly for z. 

The integral over all degrees of freedom was performed 
with the Metropolis Monte-Carlo algorithm. The ground 
state mass can be obtained most efficiently by computing 
the logarithmic derivative of G instead of G itself 



(16) 



Introducing the shorthand notation Z for the full set of 
degrees of freedom and putting S[Z] = K[z 7 s]+K[z, s] — 
V[z,z,8,s], we may write L(T) as 

L(T) = J VZ S'[Z] e- s[z] I J VZ e~ s ^ z \ (17) 

where the prime denotes an analytical differentiation of 
the functionals with respect to the endpoint T. Accord- 
ing to Eq. (17) the ground state mass is obtained by 
averaging S'[Z] over an ensemble generated by the ac- 
tion S[Z] for sufficiently large T. The FSR ground state 
wave function can readily be found by performing an 
additional integration of G in Eq. (||) over the spatial 
relative components r = y — y of the final point and 
incorporating this coordinate in the set Z. By keeping 
track of the distribution of |r|'s when computing L(T), 
the r-dependence of \I/ can be determined. 

The convergence in was studied and the mass of 
the bound state was found to become independent of N 
at typical values of N = 35-40. Furthermore, mT = 40 
usually sufficed for L(T) to become independent of T and 
to reach its asymptotic estimate (|l6|). Since the integrals 
over s and s in Eq. (j^) formally diverge for large values, 
a cutoff s max had to be introduced in order to render the 



3 



functional integrals finite. No dependence on the value 
of Smax was observed. 

In Fig. [I] we present calculations of the ground state 
mass as a function of the conventional (dimensionless) 
coupling constant g 2 /Anm 2 for the case fi/m = 0.15. 




12 3 
g /4%m 

FIG. 1. Ground state mass mo of the </? 2 x-theory as a 
function of the dimensionless coupling constant g /4irm for 
fi/m — 0.15. The inset shows the evolution of the Gross 
ground state and its unphysical branch over a larger range of 
couplings. 

Since the self-energy contributions have been neglected 
in the FSR calculations, we may directly compare the 
predictions to those of the ladder BSE and the various 
QPEs. The range of validity of the ladder theory is seen 
to be restricted to the region of small couplings. Gener- 
ally speaking, for stronger couplings all approximations 
tend to underbind the system as compared to the FSR 
results. All QPEs generate more binding energy than the 
ladder BSE and their results are generally closer to the 
FSR ones. For the Gross equation we also performed a 
calculation where the retardation in the interaction was 
neglected, i.e., we simply put u> — uj' — in the potential 
(Q) . From Fig. [I] we see that in this case the retardation 
leads to additional attraction. Particularly the ET ap- 
proximation is seen to give results that relatively provide 
the best correspondence with the FSR ones. 

We remark that due to the energy dependence in the 
two-body propagator, the Gross equation allows for a sec- 
ond, unphysical solution that starts at y/s — for g 2 = 
and for which yfs grows with increasing g 2 . This feature 
is an artefact of this particular approximation and has 
also been observed in other but similar dynamical equa- 
tions ||,^l). Inclusion of negative energy propagation 
effects was seen to cure this pathological effect. Both the 
physical and the unphysical solutions are shown in the 
inset of Fig. [j] and it is seen that they 'annihilate' each 
other at g 2 /Anm 2 ~ 5.1, for which y/s ~ 1.4m. 

In order to compare the FSR ground state wave func- 
tion to those of the ladder BSE and the various QPEs, we 



adjust the coupling constants such that the same value 
of the ground state mass is found. In Fig. || we show 
the ladder BSE and FSR wave functions for relative time 
t = and compare them to the QPE wave functions. 
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FIG. 2. Equal-time FSR and ladder BSE Euclidean wave 
function compared to those of the various QPEs. All solutions 
correspond to a bound state at mo = 1.882m. 

For convenience, the FSR wave function is normalized 
according to the standard nonrelativistic one. The mass 
of the ground state for all cases is mo = 1.882m and 
we take [i/m — 0.15. At large separations we expect 
that the wave function behavior is essentially determined 
by the binding energy of the composite system. This is 
in agreement with the calculated results shown in Fig. 
|. The main difference between the QPE predictions for 
short distances is due to the asymptotic behavior of their 
2-particle free propagator S'qpe(q) for large values of q. 
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FIG. 3. FSR and ladder BSE Euclidean wave functions at 
three values of the Euclidean relative time t. The ground 
state mass for both calculations was at mo = 1.882m. 
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Effects of relativity in the dynamics are anticipated to 
play an important role in the relative time t dependence 
of the wave function, especially at small spatial separa- 
tions between the constituents. In Fig. || we compare the 
ladder BSE and the FSR ground state wave functions for 
three values of the Euclidean relative time t. From this 
we see that the ladder BSE prediction falls off consider- 
ably faster in t as compared to the FSR result. This may 
be due to the fact that we need a substantially larger 
coupling constant in the BSE case to obtain the same 
binding energy. As a result the relativistic effects are 
enhanced in the interaction. At large distance both cal- 
culated wave functions agree essentially with each other 
and moreover show a very slow fall off in t, consistent 
with our expectation. 

For actual hadronic systems the complication of fermi- 
ons has also to be considered. Some progress has been 
achieved recently in including spin degrees of freedom 
within the FSR approach. It is clearly of great inter- 
est to study this further. In this paper we have pre- 
sented for the scalar case the first calculations of bound 
state properties beyond the ladder approximation using 
the Feynman-Schwinger representation. When compar- 
ing our results to those of the Bethe-Salpeter equation 
in the ladder approximation, we find that the crossed- 
ladders significantly contribute to the binding energy. 

It is a pleasure to thank Yu. A. Simonov for many valu- 
able and illuminating discussions concerning this work. 
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